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The modulation of charge density and spin order in (LaMn03)2n/(SrMn03)n (n=l-4) superlat- 
tices is studied via Monte Carlo simulations of the double-exchange model. G-type antiferromagnetic 
barriers in the SrMnOs regions with low charge density are found to separate ferromagnetic LaMnOs 
layers with high charge density. A metal-insulator transition perpendicular to superlattices with in- 
creasing n is observed, which provides insight into how disorder-induced localization may give rise 
to the metal-insulator transition occurring at n—3 in experiments. 

PACS numbers: 75.47.Lx, 71.30.+h, 73.21. Cd 



Introduction. Transition-metal oxide heterostructures 
provide a new avenue to utilize the complex proper- 
ties of strongly correlated electronic materials to pro- 
duce multifunctional devices. Several exotic phenomena 
emerge in these heterostructures due to the reconstruc- 
tion at the interfaces, such as the existence of a con- 
ducting state between two insulators in LaAlOa/SrTiOa 
and LaTiOa/SrTiOa. 1 As one of the most representative 
families of strongly correlated oxide materials, the man- 
ganites can also be prepared into heterostructures with 
other oxides, such as cuprates, and they exhibit interest- 
ing behavior, such as orbital reconstruction^ 

Even without involving other oxides, manganites 
heterostructures can be prepared utilizing mangan- 
ites with different doping, e.g. LaMn03 (LMO) 
and SrMn0 3 (SMO)^^i 8 These two manganites 
are parent compounds for the wide-band manganite 
Lai-^Sr^MnOs (LSMO). At low temperature (T), bulk 
LaMnOs is an A-type antiferromagnetic (A-AFM) insu- 
lator, while SrMnOs is a G-type antiferromagnetic (G- 
AFM) insulator. 9 The alloy-mixed LSMO is a ferromag- 
netic (FM) metal at low T and 0.17<x<0.5. However, 
the LMO-SMO superlattices can behave differently from 
bulk LSMO even with the same average charge density: 
(i) the ordered A-site cations in the superlattices remove 
the A-site disorder, which is important in alloy mangan- 
ites; (ii) the artificially modulated A-site cations also 
modulate the physical properties, such as charge density, 
magnetism, and conductivity. In fact, recent experiments 
on (LMO)2n/(SMO) n superlattices highlighted the ex- 
istence of an exotic metal-insulator transition (MIT) at 
n=3£&£& Moreover, LMO thin films on a SrTi0 3 (STO) 
substrate were found to be FM instead of A-AFM. 7,8 

Theoretically, in addition to ab-initio calculations, 10 
most previous model Hamiltonian investigations on 
manganite heterostructures were based on the one- 
orbital model, 11 missing the important orbital degree 
of freedom. Although more realistic two-orbital mod- 



els were used very recently ^ 12 ! 13 several properties of the 
(LMO)2n/(SMO) n superlattices are still not understood, 
particularly the explanation for the n=3 MIT. 

Models and Techniques. The two-orbital double- 
exchange (DE) model is here used to study 
(LMO) 2n /(SMO) n superlattices via Monte Carlo 
(MC) simulations. This model Hamiltonian has been 
extensively studied before and it is successful to re- 
produce the several complex phases in manganites. 14 
Details about the Hamiltonian and MC technique can 
be found in previous publications. 15 Schematically, the 
Hamiltonian reads as: 

H = H DE (t ) + #se( Jaf) + #ep(A) + - n)m, (1) 

i 

where iJoE, #se, and J/ep are the standard two- 
orbital large-Hund-coupling DE, superexchange (SE), 
and electron-phonon (EP) interactions, respectively* 1 ^ 
rii is the e g charge density at site i. \i is the uniform 
chemical potential, and is the on-site effective poten- 
tial generated by long-range Coulomb interactions that 
cannot be neglected in superlattices involving different 
electronic compositions. There are four main input pa- 
rameters: the SE coupling Jaf, the EP coupling A, /i, 
and €i. All these parameters are in units of to, which is 
the DE hopping between nearest-neighbor (NN) d 3z 2_ r 2 
orbitals along the z direction ] 14 i 15 The constant-density 
phase diagram is determined by Jaf and A. The ex- 
pected e g charge density is obtained by tuning \i. Due 
to the valence difference between La 3+ and Sr 2+ , the on- 
site Coulomb potential is inhomogeneous for the Mn 
sites. In almost all previous model investigations, the 
Coulomb interaction is treated using the Hartree-Fock 



(HF) approximation, 



11.12.13 



where the Poisson's equa- 



tion must be solved self-consistently. However, this self- 
consistent process is rather difficult to converge for the 
three-dimensional two-orbitals model when both clas- 
sical spins and lattice distortions are also MC-time evolv- 
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FIG. 1: (Color online) (a) Superlattice unit studied here, (b) 
On-site potential used in our simulation, (c) The e g charge 
density, NN spin correlation, and conductivity (see Ref.ll7lfor 
its unit) of a 4x4x4 lattice vs. /i, at T=0.01. All a are set 
to zero to simulate bulk clean- limit LSMO. 

ing. For this reason, here we adopt another strategy. 
Each €i will be determined by its eight NN A-site cation 
neighbors. 16 More specifically, in the LMO-SMO super- 
lattices, €i is for those Mn between two LaO planes 
(LMO region), it becomes V/2 for those between LaO 
and SrO planes, and finally it is V for those between two 
SrO planes (SMO region) (Fig. l(a-b)). Therefore, this 
(positive) constant V is the only parameter to regulate 
the Coulomb potential, and it is related with the dielec- 
tric constant in the HF approach. Our approximation 
is expected to capture the main physics in the LMO- 
SMO superlattices since the effective Coulomb potential 
is mainly caused by the modulation of A-site cations, and 
we believe that our qualitatively simple results shown be- 
low do not depend on these assumptions. The Coulomb 
screening by the e g electron redistribution is also taken 
into account in part by regulating the value of V. In our 
simulation, four sets of V: 0.3, 0.6, 0.9, and 1.2 are used, 
covering the realistic potential drop range between LMO 
and SMO. The important practical fact is that this ap- 
proximation enables the MC simulation on large enough 
lattices. The effort presented here would have been im- 
possible if the Poisson's equation were to be solved iter- 
atively during the MC evolution. 

For our studies, we use three-dimensional (3D, L x x 
L y x L z ) clusters with periodic boundary conditions 
(PBCs). Both L x and L y are set to 4, while L z equals 
3n (l^n^4). Thus, the superlattices are grown along the 
z direction ([001]) (Fig. l(a-b)). The MC simulation on 
the 4 x 4 x 12 lattice is already at the cutting-edge of 
current computational resources. To characterize physi- 
cal properties, the charge density, spin structure factors, 
and conductivity are calculated ] 15 ! 17 

LaMn0 3 : FM vs. A-AFM. Before the simulation of 



superlattices, it is essential to understand why the LMO 
thin films on STO are experimentally found to be ferro- 
magnetic, instead of A-AFM. For most i?Mn03 (i?=La, 
Pr, Nd, Sm, and Eu), the A-AFM phase is the ground 
stated However, in the previously obtained theoretical 
phase diagram for the two-orbitals DE model, the A- 
AFM regime was found to be rather narrow in parame- 
ter space, while the FM and orbital-ordered (OO) phase 
was clearly more robust^ Thus, to understand the FM 
nature of LMO thin films, we should consider lattice dis- 
tortions in real manganites. In the bulk, the LMO lattice 
transits from a cubic perovskite to an orthorhombic one 
at T ~ 800 K, below which the lattice constant along the 
c axis shortens compared with those along the a and b 
axes. 20 For instance, at T=300 K, l c is only ~ 0.964 Z a &, 
where l c (l a b) is the NN Mn-Mn distance along the c axis 
(within the a-b plane). Using an empirical formula by 
Zhou and Goodenough, 21 the AFM exchange intensity 
along c becomes about 1.3 times that on the a-b plane. 
Thus, this stronger AFM coupling along c will favor the 
A-AFM state. 19 However, for LMO thin films on STO, 
the LMO lattice is compressed in the a-b plane but it 
is elongated along the c axis, leading to an almost cu- 
bic crystal structure/ Therefore, the theoretical phase 
diagram, 19 derived assuming lattice isotropy, should be 
applicable to the LMO thin films that prefer the FM/OO 
phase instead of the A-AFM one. 

Following the theoretical phase diagram, 19 here we 
choose a particular set of parameters ( Jaf=0.09, A=1.2) 
for the simulation below. To justify this choice, first we 
perform a MC calculation on a 4 x 4 x 4 lattice (PBCs) 
with all €i = to examine whether the above param- 
eters are suitable for LSMO. The e g charge density n^, 
NN spin correlation (< Si • Sj >), and conductivity are 
calculated at low T (T=0.01, ~ 60 K if t ~ 0.5 eV), 
as shown in Fig. 1(c). For rii ~ 1, the phase is found 
to be FM but insulating, in agreement with the results 
found for experimental LMO thin films. The energy gap 
at the Fermi level is 0.3 (~ 150 meV), also in agreement 
with the experimental excitation energy (« 125 meV)* 7 - 
The insulating character of the state is caused by orbital 
ordering driven by Jahn- Teller distortions, although the 
orbital ordering in strained cubic LMO is a little different 
from the well-known d 3x ~ r /d 3y ~ r pattern in bulk. 22 
In the density range 0.5 < rii < 0.8, the system is FM 
and conducting, also in agreement with the LSMO prop- 
erties at the corresponding dopings. For rii < 0.5, the 
NN spin-correlations turn out to be negative, suggesting 
an AFM phase. The conductivity becomes poor with de- 
creasing rii until it reaches an insulating state. These 
rii < 0.5 behavior is also compatible with LSMO at the 
corresponding doping. Then, as a conclusion, the set 
Jaf=0.09 and A=1.2 should be a proper parameter set 
to describe (cubic) LSMO. In the following, we will use 
this set to study the (LMO)2 n /(SMO) n superlattices. 

(LMO)2n/(SMO) n superlattices. We will focus first on 
the low-T (= 0.01) MC results which can shed light on 
the ground state properties. The local e g charge densi- 
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FIG. 2: (Color online) (a) Charge modulation in the super- 
lattices with different V's (the case of L2S1 with V—1.2 can 
not been obtained due to phase separation). The two hori- 
zontal lines denote 0.5 and 2/3. (b) In-plane spin structure 
factor for V=0.9. In (a) and (b), pink bars denote SrO layers 
in the (LMO)2n/(SMO) n superlattice, while LaO layers are 
not highlighted. The lattices are simply repeated along the z 
direction if their periods are shorter than 12. 

ties of each layer vs. layer index are shown in Fig. 2(a). 
For (LMO) 2 /(SMO)i (L2S1), the charge distribution is 
fairly uniform despite the use of substantial values for 
V (0.3 — 0.9), namely the local charge density fluctuates 
weakly around the average 2/3. This is easy to under- 
stand since there is no SMO region in the L2S1 super- 
lattice and, thus, the potential fall here is only V/2. For 
other superlattices, and usually (V > 0.3), the densi- 
ties in the SMO regions are lower than 0.5. The V=1.2 
case already restricts most of the e g electrons to be in 
the LMO regions, while V = 0.3 is low enough that it 
spreads e g electrons to the SMO regions. Therefore, it 
is reasonable to conclude that the potential amplitude 
range used in our model is the proper one to cover the 
potential drop in real manganites. In the following, we 
will focus on the case V=0.9, for which the results are 
similar to the experimental data. 

Figure 2(b) shows the x-y plane spin structure factors 
vs. layer index, for F=0.9 and T=0.01. There are only 
two main components: FM and G-AFM. The latter ex- 
ists only in the SMO regions, while FM dominates in the 
LMO regions, and at the LMO/SMO interfaces. This 
spin arrangement agrees with experiments A 7 - Our calcu- 
lated spin order supports the idea that the local phases 
in superlattices are mainly determined by the local den- 
sities niocai :13 if ^locai > 0.5, the spin order is FM, and it 
is G-AFM when ni oca i~0. Therefore, other spin orders, 
e.g. A-AFM or C-AFM, may emerge once the u\ OC3 \ is 
slightly lower than 0.5. In fact, we observed the coex- 
istence of several complex spin orders at, e.g., the SMO 
region (n local ^0.15) for L8S4 with F=0.9 (Fig. 2(b)). 
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FIG. 3: (Color online) (a) In plane conductivity for the su- 
perlattices studied here, (b) Perpendicular conductivity, (c) 
Sketch of the spin order at interfaces for L4S2 (left) and L6S3 
(right), (d) Sketch of experimental setup for resistance mea- 
surements. Pink bars are SMO regions. Typical conducting 
paths via the DE process (black curves) connect NN inter- 
faces, but they will be broken when n ^ 3. 

Understanding the critical n for the MIT. One of 
the most important experimental discoveries in the 
(LMO) 2n -(SMO) n superlattices is the MIT with increas- 
ing n. To try to understand this phenomenon, here two 
conductivities (only from the DE process) are calculated 
vs. T: the in-plane one (along x or y direction) and 
the perpendicular one (along z direction), as shown in 
Figs. 3(a) and (b). All in-plane conductivities are robust 
and increase with decreasing T, suggesting metallicity. 
Our result agrees with previous studies showing that the 
charge transfer at interfaces between Mott/band insula- 
tors can generate conducting interfaces ill In contrast, 
the perpendicular conductivities show metallic behavior 
when n ^ 2, but insulating behavior when n ^ 3. To un- 
derstand this MIT, the spin arrangement at the FM/G- 
AFM interface should be considered. The intra-layer NN 
spin correlations (not shown here) confirm that spins at 
the FM/G-AFM interfaces are almost collinear at low 
T, namely the NN spins are parallel or antiparallel, as 
shown in Fig. 3(c). Therefore, when there is only one G- 
AFM layer in each superlattice unit (n=2), the "green" 
spin-up channels of the G-AFM layer link the NN FM 
layers, allowing for a good conductance. However, once 
the SrO layers thickness is 3, the two G-AFM layers cut 
off the same-spin channels, giving rise to an insulating 
behavior along the z direction. Therefore, it is natural 
to expect that n=3 must be the MIT critical point in 
(LMO) 2n /(SMO) n superlattices. 

However, it should be noted that the experimen- 
tal resistance measurements were performed using the 
four-point method (Fig. 3(d)). 6,7 Thus, to fully under- 
stand the experimental MIT, Anderson localization ef- 
fects should also be taken into consideration. When 
n ^ 3, as mentioned before, the G-AFM insulating bar- 
riers cut down conduction channels along the z direc- 
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FIG. 4: (Color online) (a) Orbital occupation of the 6th-9th 
layers for L6S3 when V = 0.9 and T = 0.01. The 6th and 9th 
are interface layers between the LaO and SrO layers, while 
the 7th and 8th are within the SMO region. Here the circle's 
area is proportional to the local e g charge density, (b) Sketch 
of orbital ordering at the interface for n—2 and 3. 

tion. Thus, the DE conducting process becomes exclu- 
sively 2D-like (i.e. only within the x — y plane). In this 
case, insulating behavior is likely induced by interface 
disorder, such as roughness, and Anderson localization, 
as already discussed in the experimental literature. 7 But 
a 3D metallicity in the n ^ 2 superlattices, avoiding the 
Anderson localization mechanism, can be achieved by 
the DE process perpendicular to the interfaces, as shown 
schematically in Fig. 3(d). In short, although the in- 
plane conductivity is metallic in our small lattices with- 
out disorder, the critical n for the MIT in the real super- 
lattices still corresponds to the same value found in our 
simulations via studies of the perpendicular conductivity. 
Orbital order. The spin configuration in Fig. 3(c) can 



also induce orbital order near the interfaces. It should be 
pointed out that the strain/stress due to the lattice mis- 
match between substrates and LMO/SMO can induce 
orbital order even without charge transfer. 10 However, 
here we propose that the the parallel spin channels (in 
both the uncut and cut cases of Fig. 3(c)) can provide 
another driving force for orbital order in manganite su- 
perlattices. Our simulations show in the n=2, 3, 4 cases, 
that the e g electrons have more tendency to occupy the 
d 3z 2_ r 2 orbital than the d x i_ y i one, to improve the ki- 
netic energy in the growth direction of the superlattice. 
One example is shown in Fig. 4. This orbital order op- 
timizes the DE process between the FM regions via the 
G-AFM layers. But in real superlattices the combined ef- 
fect of this spin-driven tendency to orbital order and that 
induced by strain/stress may compete and more quanti- 
tative calculations will be needed to decide which orbital 
order dominates. 

Conclusions. We have performed a MC simulation 
to study the two-orbitals double-exchange model for 
the (LMO)2 n -(SMO) n superlattices. First, we have ex- 
plained why the LMO thin films on STO are FM instead 
of A-AFM. Then, our simulations have shown that the 
spin order in the SMO regions is G-AFM, while it is FM 
elsewhere. The spin arrangement between the FM and 
G-AFM layers causes the metal-insulator transition, with 
n=3 as the critical value. 
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